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Abstract 



CO ■ We present CO( J = 3 — 2) emission observations with the Atacama Submillimeter Telescope Experiment 

(ASTE) toward the 5' x 5' (or 6.6 x 6.6 kpc at the distance D = 4.5 Mpc) region of the nearby barred 
spiral galaxy M 83. We successfully resolved the major structures, i.e., the nuclear starburst region, bar, 
and inner spiral arms in CO (J = 3 — 2) emission at a resolution of 22" (or 480 pc), showing a good spatial 
coincidence between CO( J = 3 — 2) and 6 cm continuum emissions. We found a global CO( J = 3 — 2) 
luminosity ^co(3-2) °^ ^-1 x 10 s K km s _1 pc 2 within the observed region. We also found ^co(3-2) m * nc 
disk region (0.5 < r < 3.5 kpc) of 4.2 x 10 8 K km s _1 pc 2 , indicating that CO(J = 3 — 2) emission in the 
disk region significantly contributes to the global ^co(3-2)- From a comparison of a CO(J = 3 — 2) data 
■ with CO (J = 1 — 0) intensities measured with Nobeyama 45-m telescope, we found that the radial profile 

of CO(J = 3 — 2)/CO(J =1 — 0) integrated intensity ratio i?3_2/i-o is almost unity in the central region 
(r < 0.25 kpc), whereas it drops to a constant value, 0.6-0.7, in the disk region. The radial profile of star 
|— j \ formation efficiencies (SFEs), determined from 6 cm radio continuum and CO( J =1 — 0) emission, shows 

the same trend as that of i?3_2/i-o- At the bar-end (r ~ 2.4 kpc), the amounts of molecular gas and the 
massive stars are enhanced when compared with other disk regions, whereas there is no excess of i?3_ 2/1-0 
and SFE in that region. This means that a simple summation of the star forming regions at the bar-end 
and the disk cannot reproduce the nuclear starburst of M 83, implying that the spatial variation of the 
dense gas fraction traced by -R3-2/1-0 governs the spatial variation of SFE in M 83. 
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1. Introduction 

In the central regions of disk galaxies, we often find in- 
tense star formation activities (e.g., Ho et al. 1997) where 
not only the star formation rates (SFRs) but also the star 
formation efficiencies (SFEs), defined as the fraction of 
the SFR surface density in the surface mass density of 
molecular g also enhanced. These star formation 

activities are referred to as starburst phenomena, which 
are different from the star formations in the disk regions 
of galaxies (Kennicutt 1998b). The starburst phenom- 
ena arc prominent not only because they show high SFRs 
but also because they show elevated SFEs. What causes 
these "high SFE" star formations? In other words, why 
do SFEs differ among galaxies (e.g., Young et al. 1996) 
or within locations/regions in a galaxy (e.g., Gao et al. 



2001)? Since stars are formed by the contraction of molec- 
ular gas, which would be initiated by gravitational insta- 
bility and/or cloud-cloud collision of molecular clouds, it 
can be easily expected that some properties of molecular 
gas play a key role in controlling the variations of SFEs. 
Investigations of molecular gas arc widely performed us- 
ing CO (J = 1 — 0) emission, which traces the total mass 
of molecular gas. In particular, extensive surveys or map- 
pings of CO( J =1 — 0) emission toward galaxies are ex- 
amined (Young et al. 1995, Kuno et al. 2000, Sorai et al. 
2000a, Nishiyama & Nakai 2001, Nishiyama et al. 2001, 
Kuno et al. 2007). However, CO(J = 1 — 0) emission si- 
multaneously traces both diffuse/cold molecular gas and 
dense/warm molecular gas, which is directly linked to star 
formation. Therefore, CO( J =1 — 0) emission is useful to 
determine the total amount of molecular gas, but it is dif- 
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ficult to compare with SFEs quantitatively. Thus, it is 
essential to study dense molecular gas, because stars are 
formed from the dense cores of molecular clouds. 

In fact, recent studies of the dense molecular medium 
in galaxies based on the observations of HCN( J = 1 — 0) 
emission, a dense molecular gas tracer (n H2 > a few xlO 4 
cm~ 3 ) due to its large dipole moment (/zhcn = 3.0 Debye, 
whereas /^co = 0.1 Debye), demonstrate the intimate con- 
nection between dense molecular gas and massive star 
formation in galaxies. For example, after the first re- 
port by Solomon et al. (1992), a tight correlation between 
HCN( J =1 — 0) and FIR continuum luminosities has been 
shown among many galaxies including nearby ones (Gao & 
Solomon 2004a, Gao & Solomon 2004b) , and high redshift 
quasar host galaxies (Carilli et al. 2005). This suggests 
that SFEs are governed by the fraction of dense molecu- 
lar gas to the total amount of molecular contents traced 
by the HCN( J= 1 -0) /CO( J= 1-0) intensity ratios (e.g., 
Solomon et al. 1992, Kohno et al. 2002, Shibatsuka et al. 
2003, Gao & Solomon 2004a). A spatial coincidence be- 
tween dense molecular gas traced by HCN( J =1 — 0) and 
massive star forming regions was also reported (Kohno et 
al. 1999, Shibatsuka et al. 2003). 

However, HCN( J = 1 — 0) emission is often too weak 
to trace a wide-area map in the disk regions of galaxies. 
HCN( J =1—0) intensity is typically 10-20 times lower 
than the CO( J =1 — 0) intensity in the disk regions of the 
Galaxy (Heifer & Blitz 1997b) and galaxies (e.g., Kohno 
et al. 1996, Heifer & Blitz 1997a, Curran et al. 2001, Sorai 
et al. 2002). Therefore, HCN( J =1-0) emission is mainly 
used to measure the amount of dense molecular gas in the 
central regions of galaxies or to obtain information inte- 
grated over the whole regions of galaxies due to inadequate 
sensitivity of current observing facilities (e.g. Nguyen- Q- 
Rieu et al. 1989, Nguycn-Q-Rieu et al. 1992, Hcnkel et 
al. 1990, Solomon et al. 1992, Israel 1992, Heifer & Blitz 
1993, Aalto et al. 1995, Curran et al. 2000, Kohno et al. 
2003). 

Instead, CO (J = 3 — 2) emission is a promising probe 
of dense molecular clouds even in the disk regions of 
galaxies. The critical density for CO(J = 3 — 2) transi- 
tion is a few x 10 4 cm" 3 . Therefore, the CO( J = 3-2) 
emission can also trace the dense component of molec- 
ular clouds as well as HCN(J = 1 — 0). In fact, in the 
central regions of nearby star forming galaxies, a tight 
correlation between CO( J = 3 — 2) and Ha luminosities, 
which is better than the correlation between CO( J — 1—0) 
and Ha luminosities, has been reported (Komugi et al. 
2007). In addition, Yao et al. (2003) discovered a trend 
that CO(J = 3 - 2)/CO(J =1-0) ratios, hereafter re- 
ferred to as i?3_2/i-0j increase with the SFE in infrared- 
luminous galaxies, which is very similar to a finding based 
on HCN( J = 1 - 0)/CO(J =1-0) ratios (Solomon et al. 
1992, Gao & Solomon 2004a). Moreover, there seem to be 
detectable amounts of CO (J = 3 — 2) emitting dense gas 
in the disk regions of galaxies such as M 51 (Wiclcbinski 
et al. 1999, Dumke et al. 2001), M 31 (Tosaki et al. 2007), 
NGC 6946 (Walsh et al. 2002), and Antennae galaxies 
(Zhu et al. 2003), where the CO(J = 3 — 2) intensities are 



found to be comparable to the CO( J =1—0) intensities. 
These results clearly indicate a potential for the use of 
CO(J = 3 — 2) emission as a dense gas tracer in the disk 
regions of galaxies. 

Spatially-resolved wide-area CO( J = 3 — 2) observations 
of nearby galaxies are also essential for the appropriate 
interpretations of high- J CO observations of high redshift 
objects (see a review by Solomon & Vanden Bout 2005), 
because most of these high- J CO observations of distant 
objects are spatially unresolved. 

Previous CO( J = 3 — 2) observations of nearby galax- 
ies, however, were almost limited to their central regions 
or to the whole regions of galaxies without resolving the 
structures using a large observing beam (e.g. Wild et al. 
1992, Hurt et al. 1993, Devereux et al. 1994, Mauersberger 
et al. 1996, Wielebinski et al. 1999, Nieten et al. 1999, 
Mauersberger et al. 1999, Harrison et al. 1999, Sandqvist 
1999, Dumke et al. 2001, Liszt 2001, Meier et al. 2001, 
Wiedner et al. 2002, Vila-vilaro et al. 2003, Yao et al. 

2003, Hafok & Stutzki 2003, Matsushita et al. 2004, Iono 
et al. 2004, Narayanan et al. 2005, Weifi et al. 2005, Israel 
2005, Israel et al. 2006). 

In this paper, we present wide-area observations 
of CO(<7 = 3 — 2) line emission using the Atacama 
Submillimctcr Telescope Experiment (ASTE; Ezawaet al. 

2004, Kohno 2005) towards the whole inner disk of M 83 
(NGC 5236). M 83 is a nearby, face-on, barred, grand- 
design spiral galaxy hosting an intense starburst at its 
center. The distance to M 83 is estimated to be 4.5 Mpc 
(Thim et al. 2003); therefore, 1" corresponds to 22 pc. 
The inclination of M 83 is 24° (Comte 1981). Its proxim- 
ity and face-on view enable us to resolve its major struc- 
tures such as the nuclear starburst region, bar, and inner 
spiral arms even by singlc-dish observations at millimeter 
or submillimctcr wavelength. For these reasons, M 83 is 
the best target to investigate the spatial variations of SFE 
and the physical state of the molecular gas between the 
nuclear region and the disk region. The CO observations 
toward M 83 were reported for various transitions (see a 
comprehensive summary in Lundgrcn ct al. 2004) . Handa 
et al. (1990) present a CO( J =1 — 0) map of the nucleus 
and the bar at a resolution of 16" using a Nobcyama 45- 
m telescope. A new wider CO( J =1 — 0) map using the 
Nobeyama 45-m telescope is now available (Kuno et al. 
2007). Moreover, Crosthwaite et al. (2002) and Lundgren 
et al. (2004) present CO(J= 1-0) and CO(J = 2-l) maps 
covering the entire optical disk. However, the high- J CO 
line observations such as CO (J = 3 — 2) were still limited 
to the nuclear region and the bar (Petitpas & Wilson 1998, 
Israel & Baas 2001, Dumke et al. 2001, Sakamoto ct al. 
2004, Bayet et al. 2006). 

The goals of this paper are (1) to depict the global dis- 
tribution of CO( J = 3 — 2) line emission or dense molec- 
ular gas in the disk region of M 83, (2) to investigate the 
spatial variations of i?3_2/i-o within the disk of M 83, 
by computing the radial profiles of CO( J = 3 — 2) and 
CO (J =1 — 0) intensities, and (3) to examine the corre- 
lation of the A3-2/1-0 values and SFEs through a com- 
parison of the radial profiles. In particular, we focus on 
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Right ascension (J2000.0) 

Fig. 1. Points observed (crosses) using the ASTE super- 
posed on a V-band image of M 83 obtained with VLT 
(Comeron 2001). The grid points are spaced by 11" at the 
center, bar, and along the spiral arms and by 22" in the in- 
ter-arm. 

the difference of properties of the nuclear starburst and 
star formation in the disk region. Further analysis of our 
CO( J = 3 — 2) data of M 83 will be presented in forthcom- 
ing papers. 

2. Observations and Data Reduction 

CO( J = 3 — 2) observations towards M 83 were per- 
formed using the ASTE from December 2004 to January 
2005, and June 2005. The size of the CO( J = 3-2) map 
is about 5' x 5' (6.6 x 6.6 kpc), including the nuclear star- 
burst region, the bar, and the inner spiral arms. We set 
the grid spacing to 11", except for the inter-arm region 
where the grid spacing was set to 22". We set the posi- 
tion angle of the major axis to 45° (Comte 1981). The 
total number of points observed is 419, and the observed 
positions are indicated in figure 1. The full-half power 
beam width (HPBW) is 22", which corresponds to about 
480 pc at a distance of 4.5 Mpc. 

We used a cartridge-type cooled SIS mixer receiver, 
SC345, developed by collaboration among the University 
of Tokyo, the Osaka Pref. University, the Nagoya 
University, and the National Astronomical Observatory 
Japan (NAOJ). The receiver SC345 was plugged into the 
Single Cartridge Dewar of the ASTE, which is based on 
the design of the ALMA test cryostat (Sekimoto et al. 
2003, Sugimoto ct al. 2003, Yokogawa et al. 2003) using a 
3-stage Gifford McMahon cryocooler with an He pot for 
temperature stabilization. A Nb-based parallel-connected 
twin-junction-type mixer is mounted on the 4 K stage of 
the dewar with a corrugated feed horn and a local oscilla- 
tor (LO) signal coupler to inject radio-frequency (RF) and 
LO signals into the receiver. The intermediate-frequency 
(IF) range is 4.5 to 7.0 GHz. An LO system is com- 
posed of a chain of multipliers, i.e., an active frequency 



doublcr, two frequency doublers, and a frequency triplcr, 
provided by Virginia Diode Inc. A source signal generated 
by a GPS-locked frequency synthesizer with a frequency 
range of 13.5 GHz to 15.5 GHz is multiplied by 24 us- 
ing the multiplier chain system. Both the attenuation of 
the LO signal power and the change in the frequency set- 
ting without a phase lock loop of Gunn oscillator easily 
facilitate fully remote operations of the receiver system. 
The resultant LO frequency range is 324 to 372 GHz, and 
the typical receiver noise temperature is about 100 K (in 
double side-band) for this frequency range. The average 
system temperature was 200 K (in double side-band) dur- 
ing our observations of M 83. The temperature scale was 
calibrated by the chopper-wheel method, using a room 
temperature single absorber in front of the receiver. The 
observations were made remotely from the ASTE oper- 
ation rooms of the Institute of Astronomy (IoA), Univ. 
of Tokyo, NAOJ, and the Nobeyama Radio Observatory 
(NRO) using a network observation system N-COSMOS3 
developed at NAOJ (Kamazaki et al. 2005). The sky emis- 
sion was subtracted by position switching with two off- 
source positions at offsets in an azimuth of ± 10' from the 
central position, which was taken at a = 13 h 37 m 00 s .48, 
5 = -29°51'56".5 (J2000). This position corresponds to 
the radio continuum center (Sukumar et al. 1987). The 
spectra of the CO( J = 3 — 2) emission were obtained with 
four units of XF-type digital spectrometers, MAC (Sorai 
et al. 2000b), with a band width of 512 MHz and 1024 
channels, corresponding to a velocity coverage of 440 km 
s _1 for each unit. The absolute pointing of the antenna 
was checked every hour using Jupiter and the accuracy 
was found to be better than ~ 2" rms. 

We observed the CO( J = 3-2) emission of IRC+10216 
twice a day to verify the stability of the main beam 
efficiency (tjmb) of the ASTE 10-m dish, and we have 
compared our data with the CO( J = 3 — 2) emission of 
IRC+10216 obtained by CSO observations (Wang et al. 
1994), which are carried out using an SSB-filter. For the 
CSO observations, the line peak temperature of CO( J = 
3-2) emission of IRC+10216 in T MB was 32.5 ± 0.5 K. 
Our CO(J = 3 — 2) data was calibrated using this value. 
The resultant ?7mb obtained by the calibration of the line 
temperature was 0.63 ± 0.06. In addition to this, we have 
also calculated tjmb using the brightness temperature of 
Jupiter recorded during the pointing observation. We as- 
sumed that the brightness temperature of Jupiter at 345 
GHz is 174 K (Mangum 1993). The resultant t^mb de- 
duced from the Jupiter measurements was 0.62 ± 0.03. 
These ?7mb values obtained by the two methods are con- 
sistent, indicating that the side-band gain ratio of the re- 
ceiver SC345 was almost 1. 

Data reduction were carried out using NEWSTAR and 
GBASE, tools for single-dish data analysis developed by 
NAOJ. The linear baselines were removed from most of 
the spectra. We binned the adjacent channels to a velocity 
resolution of 5 km s" 1 for the CO( J = 3 — 2) spectra. The 
resultant rms noise level was typically in the range of 40 to 
70 mK in the Tmb scale. The observation parameters for 
M 83 are listed in table 1. Also, the adopted parameters 
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Tabic 1. Observation parameters of M 83. 



Telescope, Receiver 
Spectrometer 
Velocity coverage 
Map description 
rms noise level 
Beam size 

Main-beam efficiency* 



ASTE 10-m, SC345 

512 MHz, 1024 channel 

440 km s" 1 

419 positions, 5' x 5' 

40-70 mK (in T MB scale) 

22" 

0.62 ± 0.06 ± 0.03 



For the main-beam efficiency, the first error indicates 
systematic error and the second, random error. 

Table 2. Adopted parameters of M 83. 



Map center": 

R.A. (J2000) 

Dccl. (J2000) 
Distance^ 
Linear scale 

Systemic velocity (LSR) C 
Position angle d 
Inclination' 1 



13 h 37 m 00 s .48 

-29°51'56".5 

4.5 Mpc 

22 pc arcsec - 

510 km s -1 

45° 

24° 



Reference - 



"Sukumar et al. (1987), 6 Thim et al. 
c this paper, d Comte (1981) 



(2003), 



of M 83 are listed in table 2. 
3. Results 

3.1. CO (J = 3-2) Spectra 

Figure 2 shows all 419 spectra of the CO( J = 3 — 2) emis- 
sion of M 83. The X and Y axes were set parallel to the 
major and minor axes of the galaxy, respectively. In ad- 
dition, figure 3 shows the spectra in the central 44" x 44" 
region. Significant CO(J = 3 — 2) emission for which the 
signal to noise ratio exceeds 4 is detected for 275 points out 
of the 419 observed points. Strong emission is observed at 
the nucleus and the leading edges of the bar. Moreover, 
the emission line is significantly detected from most of the 
points along the inner spiral arms. No significant emission 
lines are detected in the inter-arm regions. 

The physical properties of the observed CO(J = 3 — 2) 
emission line of M 83 are summarized in table 3. Velocity- 
integrated intensities, ^co(3-2) > were calculated from the 
CO emission within a velocity range of 380 to 620 km s _1 
of spectra. In addition, we calculated the luminosity of 

was 



L' 



We found that L' r 



CO(J = 3-2), ^ CO (3-2)- vvc; "J-"*" ^CO(3-2) 

9.3 x 10 7 K km s -1 pc 2 in the central region (r < 0.5 kpc) 
and 4.2 x 10 8 K km s -1 pc 2 in the disk region (0.5 < r < 
3.5 kpc), respectively. The resultant global CO luminosity 
was found to be 5.1 x 10 s K km s -1 pc 2 . This means that 
■^CO(3-2) m ^ ne disk region significantly contributes to the 
global ^co(3-2)- Note that the global £' CO ( 3 _2) * s s * m 
the lower limit because we did not observe CO( J = 3 — 2) 
emission toward the outer disk region (r > 3.5 kpc). 



Now, we compare L 



CO(3-2) 



of M 83 to those of other 



galaxies. Yao et al. (2003) gave the plot of L' CO ( 3 _ 



2) 



Lpm for 54 galaxies of which Lfir is in the range of 10 9 Lq 
to 1O 12 L . L Fm of M 83 is ~ 9 x 1O 9 L (Mauersberger 
et al. 1999), and the average ^co(3-2) m Yao's sample is 
- 10 4 - 5 K km s" 1 Mpc 2 n b at L F ir = 9 x 1O 9 L . Because 
10 4 - 5 K km s" 1 Mpc 2 n b in Yao et al. (2003) corresponds 
to 2 x 10 8 K km s' 1 pc 2 , £' CO (3_ 2 ) of M 83 determined by 
our wide-area CO( J = 3 — 2) observations is comparable 
to that of Yao's sample. 

The peak temperature and the integrated intensity at 
the nucleus (r < 250 pc) are 1.51 ± 0.04 K and 161 ± 
4 K km s -1 , respectively, in the Tmb scale. This inte- 
grated central CO value is in excellent agreement with 
167 ± 15 K km s -1 for 21" beam reported by Israel & 
Baas (2001). Note that Mauersberger et al. (1999) report 
the CO (J = 3 — 2) integrated intensity at the nucleus of 
M 83 of 126 K km s _1 , whereas Dumke et al. (2001) found 
that the line intensity was 234 ± 14 K km s -1 , although 
they used the same 10-m telescope, the Heinrich-Hertz- 
Telescope. Our measured CO( J = 3 — 2) integrated in- 
tensity lies between these two values. It is unclear what 
causes such a significant discrepancy. There is slight dif- 
ference in the position of the "center" , a few arcsec, yet it 
is clearly insufficient to account for the large discrepancy 

At the bar-end (110" or 2.4 kpc away from the center 
to the X direction of M 83; Kuno et al. 2007), where sec- 
ondary intensity peaks exist (see also figure 4 and figure 
7a), the peak temperatures are in the range of 0.5 K to 0.9 
K. The typical line width is about 50 km s _1 . In the spiral 
arm region, the line temperature is in the range of 0.3 to 
0.5 K. In the spectra of M 51 (Wielebinski et al. 1999), 
the peak temperatures of the CO(J = 3 — 2) emission in 
the spiral arms are 0.3-0.4 K, which are comparable to 
those of M 83. 

3.2. Global Distributions of the CO (J = 3 — 2J Emission 

Figure 4 shows an integrated intensity map and a peak 
temperature map of the CO( J = 3 — 2) emission of M 83. 
Extended CO (J = 3 — 2) emission along the bar can 
be clearly seen as well as a prominent concentration of 
CO( J = 3 — 2) emission toward the nucleus. The extended 
components are in the leading sides of the bar; they are 
often referred to as "molecular offset ridges," which are 
visible in high resolution CO (J =1 — 0) maps of barred 
galaxies. A strong CO( J = 3 — 2) emission is also evident 
along the spiral arms; in particular, we can clearly trace 
two "CO( J = 3 — 2) arms" in the peak temperature map. 
Note that they are not so clear in the integrated intensity 
map due to the narrow line width and weaker intensity of 
CO in the spiral arms. 

Previous CO( J = 3 — 2) observations of M 83 were con- 
fined to a small portion (~ V — 2') of the galaxy (Pctitpas 
& Wilson 1998, Israel & Baas 2001). Dumke et al. (2001) 
mapped a wider region (~ 2' x 2'); however, detected 
CO( J = 3 — 2) emission remained confined within a re- 
gion along the bar. Thus, this is the first map that clearly 
depicts the distribution of warm and dense molecular gas 
clouds in the disk region of M 83 including spiral arms. 
Basically, the overall appearance of our CO (J = 3 — 2) 
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Fig. 2. All measured spectra of CO( J = 3 — 2) emission of M 83. The scale of the spectra is indicated by the small box inserted 
in the lower right corner. Significant emission is detected not only at the nucleus but also at the bar and inner spiral arms. The 
central 5x5 spectra are magnified and shown in figure 3. 



Table 3. Physical properties of observed CO( J = 3 — 2) emission line of M 83. 



Region 


^CO(3-2) 


^CO(3-2) 


•R3-2/1-0 




[K km s- 1 ] 


[K km s^ 1 pc 2 ] 




M 83 center (r < 0.5 kpc) 


120 


9.3 x 10 7 


0.97 ± 0.10 


disk (0.5 < r < 3.5 kpc) 


10 


4.2 x 10 8 


0.65 ± 0.07 


global 


12 


5.1 x 10 s 


0.69 ± 0.07 



G 
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Fig. 3. The central 5x5 spectra of M 83. Grid spacing is 
11". 

maps is similar to those of previous CO(J = 1 — 0) and 
CO(J = 2 — 1) maps presented by Handa et al. (1990), 
Crosthwaite et al. (2002), Lundgren et al. (2004), Kuno 
et al. (2007). Our high-quality data set allows us to 
see the molecular gas structures more clearly through 
CO( J = 3 — 2) emission. 

In figure 5, we see a striking similarity between the 
contour map of the CO(J = 3 — 2) peak temperature 
and a 6 cm radio continuum map obtained with the 
VLA (Ondrechen 1985, Neininger et al. 1993), which 
is convolved to 22" resolution to compare it with our 
CO( J = 3 — 2) image. The centimeter continuum emis- 
sion is mainly dominated by nonthermal synchrotron radi- 
ation from supernova remnants in M 83 (Ondrechen 1985, 
Turner & Ho 1994). The observed radio intensity is pro- 
portional to a recent star formation rate. We find that 
these maps are well correlated, indicating a good spa- 
tial correlation between star formation and the dense and 
warm gas traced by the CO( J = 3 — 2) emission. 

Then, we examine a comparison of CO (J = 3 — 2) in- 
tensity with CO( J =1 — 0) intensity. We use new CO( J = 
1—0) data obtained using the Nobeyama 45-m telescope 
(Kuno et al. 2007). The HPBW of the CO(J = 1-0) 
observations with the 45-m telescope is ~ 16", and the 
central CO( J =1 — 0) intensity is 170 K km s _1 in the 
Tmb scale. Here, we compare the new CO( J = 1 — 0) data 
with previous CO( J =1 — 0) data presented by Handa et 
al. (1990). According to Handa et al. (1990), the central 
CO( J =1 — 0) intensity at a resolution of 16" was ~ 80 K 
km s _1 in the T£ scale. Since the beam efficiency of their 
observations was 0.45, the resultant central CO( J =1 — 0) 
intensity in the Tmb scale was 180 K km s _1 . Therefore, 
the new CO( J =1 — 0) data is in excellent agreement with 
the previous CO( J =1 — 0) data. 

Figure 6 displays a plot of CO (J = 1 — 0) integrated 
intensity vs. CO( J = 3 — 2) integrated intensity for each 



observed point. To match the beam size of CO( J =1 — 0) 
to that of CO( J = 3-2), we convolved the CO( J = 1- 0) 
data of 16" resolution to 22" resolution. The red circles 
indicate the observed points in the central region (r < 0.5 
kpc) of M 83 and the black ones, the inner disk region 
(0.5 < r < 3.5 kpc). This plot excludes the data points 
where the CO(J= 1 — 0) intensities are lower than 8 K km 
s _1 , corresponding to a typical 3 a rms. At the center of 
M 83, the CO(J = 3 — 2) intensity is close to or greater than 
the CO(J =1-0) intensity, whereas the CO(J = 3-2) 
intensities are weaker than the CO( J =1 — 0) intensity for 
most of the data points in the disk region. 

The azimuthally averaged radial profiles of CO(J = 
3 — 2) and CO (J =1 — 0) integrated intensities are shown 
in figure 7a. The radial profile of CO( J = 3 — 2) is similar 
to that of CO (J =1 — 0); both of them have strong peaks 
at the nucleus and a secondary peak at the bar-end (~ 2.4 
kpc), where the molecular gas is accumulated. Such a ra- 
dial CO profile is often observed in barred spiral galaxies 
(Nishiyama et al. 2001). Note that the radial profile of 
CO( J = 3 — 2) integrated intensities for larger radii (1 < 
r < 3.5 kpc) are biased toward the intensites in the bar 
and the inner spiral arms since the observed points of the 
CO(J = 3 — 2) emission are biased toward the bar and 
the spiral arms. The same bias exists in the radial pro- 
file of i?3_2/i-0; which will be described in the following 
subsection. 

3.3. CO(J = 3 - 2)/CO(J = 1-0,) Intensity Ratio 
-R3-2/1-0 

We have obtained -R3-2/1-0 for each observed point. 
The absolute error of the CO(J = 3-2) and CO(J = 
1 — 0) intensities are estimated to about 15 %, respectively. 
Thus, the absolute error of is estimated to be 

about 30%. 

The central i?3_2/i-o a t a resolution of 22" is 1.1. For 
the center of M 83, Dumke et al. (2001) reported -R3-2/1-0 
at a resolution of 22" is ~ 1.4. This discrepancy is mainly 
due to their high CO (J = 3 — 2) integrated intensity value 
than our CO( J = 3 — 2) data (see section 3.1). 

Note that Kramer et al. (2005) reported that the central 
J?3_2/i_o is 0.6 at a resolution of 80". Unfortunately, 
we cannot convolved our CO (J = 3 — 2) data to 80" due 
to imperfect sampling in the disk region. Therefore, we 
cannot compare i?3_2/i-o a t a resolution of 80" with that 
of Kramer et al. (2005). 

A radial profile of -R3-2/1-0 i s displayed in figure 7b. 
i?3_2/i-o is almost unity in the nuclear region (r < 0.25 
kpc), whereas it drops to a constant value in the regions 
with larger radii (1 < r < 3.5 kpc). The mean value of 
the ratio is 0.65 ± 0.07 within the inner disk region. 

There is a discrepancy between this radial profile of 
the # 3 _2/i-o a nd the CO( J = 2 - l)/CO(J =1-0) inte- 
grated intensity ratio -R2-1/1-0 presented by Lundgren et 
al. (2004). Their -R2-1/1-0 radial profile has a peak not 
at the nucleus but at the radius of 1 kpc, corresponding 
to the bar region. The reason for this discrepancy is un- 
clear, although the difference in the observing beams (our 
i? 3 _2/i-o is obtained at a 22" resolution, whereas their 
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Fig. 4. (left) A contour map of the CO( J = 3 — 2) integrated intensity of M 83. The central cross indicates the central reference 
position of the map. Contours are at 5, 10, 15, 20, 25, 30, 40, 50, 60, 80, 100, 120, 140, and 160 Kkrns" 1 , and the peak is 161.2 
Kkms -1 . The 1 a noise level of the map is 2 Kkms -1 . (right) A contour map of the CO(J = 3 — 2) peak brightness temperature. 
Contours are at 0.2, 0.3, 0.4, 0.5, 0.7, 0.9, 1.1, 1.3, and 1.5 K, and the peak is 1.53 K. Major structures, i.e., the nuclear starburst 
region, bar, and inner spiral arms are clearly resolved. 
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Fig. 5. The CO( J = 3 — 2) peak temperature map (white 
contour) superposed on the 6 cm radio continuum map (grey 
scale) obtained with the VLA (Ondrechen 1985). These two 
maps correlate well spatially. 



i? 2 -i/i-o is & t a 49" resolution) may partly account for 
this. 

In other galaxies, the trend of the spatial variations of 
i? 2 -i/i-o is similar to that of -R3-2/1-0 hi M 83. For 
example, the reported -R2-1/1-0 m the central 900 pc of 
the Galaxy is 0.96 (Oka ct al. 1998, Sawada et al. 2001), 
which is higher than the typical value of i? 2 -i/i-o m the 
Galactic disk, 0.6-0.7 (Sakamoto et al. 1995, Sakamoto et 
al. 1997). Further, for the Large Magellanic Cloud, Sorai 
et al. (2001) reported that there existed a radial gradient 
of i?2-i/i-o °f 0-94 in the inner region (<2 kpc from the 
kinematic center) and 0.69 in the outer region. 

We divided the -R3-2/1-0 values of M 83 into two 
groups, i.e., the central region (r < 0.5 kpc) and the 
disk region (0.5 < r < 3.5 kpc). The intensity-weighted 
(summation of the intensity in each i? 3 _ 2 /i_o range) his- 
tograms of -R3-2/1-0 fo r these two groups are shown in 
figure 8. These histograms include only the data for which 
signal to noise ratio of the CO (J = 3 — 2) emission exceeds 
5. The number of observed points used to make this his- 
togram is 11 for the center and 219 for the disk. 

In the central region (r < 0.5 kpc) of M 83, ^£3-2/1-0 
is in the range of 0.7 to 1.2, and the average value is 0.97. 
This is comparable with the average value of i? 3 _ 2/1-0 in 
nearby starburst galaxies, ~ 0.9 (Devereux et al. 1994), 
and slightly higher than the average values of i?3_ 2 /i_o m 
compact galaxies, ~ 0.7 (Israel 2005); infrared luminous 
galaxies, ~ 0.6 (Yao et al. 2003); dwarf starburst galaxies, 
~ 0.6 (Meier et al. 2001); and early type galaxies, ~ 0.4 
(Vila-vilaro et al. 2003). This could be partly due to the 
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Fig. 6. A plot of the CO(J = 1-0) intensity vs. 
CO( J = 3 — 2) intensity for each observed point of M 83. Data 
points where CO( J = 1 — 0) intensities are lower than 8 K km 
s -1 are excluded. The red circles indicate the points in the 
central region (r < 0.5 kpc) and the black ones, the inner disk 
region (0.5 < r < 3.5 kpc). In this plot, the error originated 
from the noise is only included, but errors from uncertainties 
of baselines and calibrations arc not included. The red line 
represents the average i?3_ 2 /i_o in the central region, 0.97, 
while the black line represents that in the disk region, 0.65. 
At the center of M 83, the CO(J = 3 — 2) intensities are close 
to or greater than CO( J = 1 — 0), whereas the CO(J = 3 — 2) 
intensities are weaker than CO( J = 1 — 0) for most of the data 
points in the disk region. 



difference in the aperture sizes where the R 



3-2/1-0 



val- 



ues were measured (e.g., Meier ct al. 2001). Moreover, 
our i?3_ 2 /i_o in the central region of M 83, 0.97, ap- 
pears to be slightly higher than the average values of 
i?3_2/i-o m nearby galaxies, ~ 0.6 (Mauersberger et al. 
1999). However, as we pointed out in section 3.1, /co(3-2) 
at the nucleus of M 83 reported by Mauersberger et al. 
(1999) is about 30% lower than our result. Therefore, we 
cannot simply compare these two average values. The 
reported -R2-1/1-0 ar e also close to unity (~ 0.7 to 1; 
Crosthwaite et al. 2002, Lundgren et al. 2004), although 
their aperture sizes are larger than ours by a factor of 2 
or more. 

On the other hand, the ratio in the inner disk spreads 
over a wide range of 0.2 to 1.4, and the average value is 
0.65. More than half of the data points in the disk re- 
gion show i?3_2/i_o values ranging from 0.5 to 0.8. These 
i?3_2/i-o values are slightly higher than the average value 
of cold giant molecular clouds (GMCs) in the Galactic 
disk, 0.4 (Sanders et al. 1993). Note that there exist four 
data points for which -R3-2/1-0 exceeds 1.2 even in the 
disk region. These four points are in the spiral arm re- 
gions but are not adjacent to each other, i.e. existing 
locally. A similar situation has been observed in a spiral 
arm of the Milky Way galaxy. i?3_ 2 /i_o in the Sagittarius 
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Fig. 7. (a) CO( J = 3-2) and CO( J = 1-0) line intensi- 
ties as a function of the galactoccntric radius of M 83. Both 
profiles show strong peaks in the central region as well as sec- 
ondary peaks at the bar-end (r ~ 2.4 kpc). (b) ^3-2/1-0 as 
a function of the galactocentric radius of M 83. R3-2/1-0 
is almost unity at the nucleus (r < 0.25 kpc) and drops to 
a constant value, 0.6—0.7, in the disk region (0.5 < r < 3.5 
kpc), showing no secondary peak at the bar-end. 

arm of the Galaxy is ~ 0.5, but there exist locations at 
which i? 3 _2/i_o exceeds unity (T., Sawada. private com- 
munication). 

4. Discussion 

Thanks to our high-quality data set of CO(J = 3 — 2) 
emission, we have obtained the distribution of wide- 
ranging i? 3 _2/i_o, which can resolve the 500 pc scale 
structure. As mentioned in section 3, i?3_2/i-o is almost 
unity at the nucleus, and then drops to a constant value 
(0.6-0.7) in the disk regions. 

In this section, we discuss the relation between the phys- 
ical state of molecular gas and star-formation activity us- 
ing the data on i?3_2/i-o- 

4-.1. Interpretation of R 3 _2/ 1-0 

To obtain the physical interpretation of the observed 
i?3_2/i-o values, we have examined the dependence of the 
-R3-2/1-0 values on the kinetic temperature and gas den- 
sity by employing the large velocity gradient (LVG) ap- 
proximation (Scoville & Solomon 1974, Goldreich & Kwan 
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Fig. 8. Intensity- weighted histograms of Rz—2/1—0 of M 83. 
(top) at the center (r < 0.5 kpc). (bottom) in the disk region 
(0.5 < r < 3.5 kpc). 

1974). Figure 9a shows the result of the LVG calculations, 
assuming Z ( 12 CO) /dv/dr = 8x 10~ 5 , which would be a 
typical value for galaxies. If we consider the one-zone 
model to be valid and assume the kinetic temperature 
Tkin ~ 20 K, which seems to be relevant for GMCs in the 
disk regions, it appears that -R3-2/1-0 depends mainly on 
the density of gas rather than the kinetic temperature in 
the disk region where the -R3-2/1-0 is in the range of 0.1 
to 0.7. The H2 gas density is ~ 10 25-3 cm~ 3 in the disk 
region. However, in this assumption of Z( 12 CO)/dv/dr = 
8 x 10 -5 , we cannot reproduce the high i? 3 _ 2 /i_o values 
observed in the central region, where i? 3 _ 2 /i_o is close to 
or exceeds unity, suggesting a low opacity of CO emis- 
sion or difference in emitting regions of a clump between 
CO(J = 3-2) and CO( J =1-0) due to external heat- 
ing (Gierens et al. 1992). Therefore, in the central region, 
we must assume another condition different from that in 
the disk region. Figure 9b shows the result of the LVG 
calculations, assuming Z( 12 CO) /dv/dr = 8 x 10~ 6 . In this 
condition, the kinetic temperature is higher than 60 K and 
the H2 gas density is ~ 1 x 10 4 cm" 3 at the nucleus where 
the i? 3 _2/i_o exceeds unity. This is a fairly strong con- 
straint, and the molecular gas is in a significantly warm 
and dense state. From these results of the LVG approxi- 
mation, the variation of i?3_2/i_o in the range of 0.4 to 1.0 
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Fig. 9. Curves of constant R 



'3-2/1-0 



as functions of kinetic 



temperature and density of molecular gas. (a) CO fractional 
abundance per unit velocity gradient Z( 12 CO) / (dv / dr) was 
set to 8 x 10~ 5 , where Z( 12 CO) is defined as [CO]/[H 2 ], and 
the unit of dv/dr is km s -1 pc" 1 . Note that no contours 
showing iJ3_2/i_o of 1 are seen, (b) CO) / (dv / dr) was 
set to 8 X 10 -6 . 

seems to depend mainly on the density of the molecular 
gas. Note that in the condition where i?3_ 2 /i_o is close 
to unity, i?3_2/i-o also restricts the kinetic temperature 
as well as the density of the molecular gas. 

In the case of external galaxies, the observing beams are 
usually too large to distinguish the individual molecular 
cloud structures, and the one-zone assumption is no longer 
valid. Thus, i?3_2/i-o could be a rough measure of the 
fraction of the dense molecular gas (nn 2 > a f ew xlO 3 
cm~ 3 ) to the total (including diffuse) molecular clouds 
within the observing beam. 

4-2. Relation between the i?3_2/i-o an d SFE 

The observed radial variations of -R3-2/1-0 an d the LVG 
calculation suggest the fraction of the dense molecular gas 
at the nucleus is higher than that in the disk region in 
M 83. Here we investigate the physical relationship be- 
tween the molecular gas properties traced by -R3-2/1-0 
and the massive star formation. The first step is to quan- 
tify the star formation rates (SFRs) and star formation 
efficiencies (SFEs) in M 83. 
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4.2.1. Derivation of SFR and SFE 

An SFE is expressed using an SFR or a surface density 
of SFR (Ssfr) and a surface mass density of molecular 
hydrogen (Sh 2 ) as follows: 



SFE 



yr" 



J SFR 



M yr !pc 2 



(1) 



To calculate Ssfr, we used a 6 cm radio continuum 
map obtained with the VLA (Ondrechen 1985), as shown 
in figure 5. We derived a radial distribution of the 6 cm 
radio luminosity. Condon (1992) relates the nonthermal 
luminosity Ln and SFR as follows: 



SFR (M > 5M G 
A/oyr- 1 



1.9 x 10" 



( U 



VWHz- 



GHz 



(2) 



where a ~ 0.8 is a nonthermal spectral index and v is the 
observed frequency. By assuming that all the observed 6 
cm radio continuum flux of M 83 originate from nonther- 
mal emission due to supernova remnants, we converted 
the radial distribution of 6 cm radio luminosity into that 
of Ssfr by employing equation (2). This is shown in 
figure 10a. It is similar to those of CO( J = 3 — 2) and 
CO(J= 1 — 0) integrated intensities; there is a strong peak 
at the center. In addition, we can see a secondary peak 
at the bar-end (r ~ 2.4 kpc). The azimuthally averaged 
SFRs are about 1 x 10~ 7 M0yr _1 pc~ 2 in the disk region 
and about 3 x 10~ 6 M yr _1 pc~ 2 at the nucleus of M 83, 
Ssfr at the nucleus where a moderate starburst oc- 



l.e 



curs is enhanced by a factor of ~ 30 as compared with 
that in the disk region. 

Sh 2 is estimated from the CO( J = 1 
a A^h 2 / Ico conversion factor Xco- 



- 0) intensity using 
Here, we adopted 



a Galactic X C o value of 1.8 x 10 20 cm" 2 (K km s" 1 )" 1 
(Dame et al. 2001). Note that the effect of the constant 
Xqo will be discussed at the end of this subsection. £h 2 
is calculated as 



M oP c- 



= 2.89 x cos(i) 



CO 



K km s 1 



1.8 x 10 20 cm- 2 (Kkms- 1 )- 



(3) 



where Ico is the CO(J =1 — 0) intensity and i is the 
inclination of the galaxy (i = 24° for M 83, Comte 1981). 

Figure 10b shows the resultant azimuthally averaged 
SFEs as a function of the galactocentric radius using equa- 
tion (1). In contrast with the radial profiles of CO and 
Esfr, the radial profile of SFE in figure 10b shows a re- 
markable difference; there is no enhancement of SFE at 
the bar-end (r ~ 2.4 kpc). In the nuclear region, SFE is 
8.9 x 10~ 9 yr" 1 , whereas the average SFE is 2.2 x 10 -9 



yr 



in the disk region, i.e., four times lower than that 
in the nuclear region. Considering the fact that these de- 
rived SFR and SFE values are limited to contributions 
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Fig. 10. (a) The azimuthally averaged SFR as a function of 
the galactocentric radius of M 83. The profile of SFRs show 
strong peaks in the central region and secondary peaks at the 
bar-end (r ~ 2.4 kpc) as well as CO intensities shown in fig- 
ure 6a. (b) The azimuthally averaged SFE and ^3—2/1-0 as 
a function of the galactocentric radius of M 83. Both pro- 
files show a strong peak at the center, whereas no significant 
enhancement is visible at the bar-end, despite the fact that 
there is a secondary peak in the CO intensities and SFRs. 

from massive stars heavier than 5 Mq, the total (i.e., in- 
cluding low-mass stars) SFE at the nucleus would be of 
an order of a few x 10~ 8 yr -1 (see Condon et al. 2002), 
which is in fact within the range of circumnuclear star- 
bursts (Kcnnicutt 1998a). The SFEs in the disk region of 
M 83 are comparable to or slightly higher than those of 
nearby gas-rich spiral galaxies (e.g., Wong & Blitz 2002). 

Note that we adopt a constant Xqo value over the ob- 
served region, but our conclusion that there is an excess 
of SFE at the center of M 83 is not weakened even after 
considering a possible variation of Xco within M 83. In 
the central regions of various galaxies, Xco factors tend 
to show smaller values than those in the disk region by 
a factor of 2-3 or more (e.g., Nakai & Kuno 1995, Regan 
2000), including the Galactic Center (e.g., Oka et al. 1998, 
Dahmen et al. 1998). A smaller Xco yields a smaller gas 
mass in the central region, providing higher SFE than the 
present one in figure 10b. 
4-2.2. Comparison of R3-2/ 1-0 with SFE 

Now we can compare the radial profile of -R3-2/1-0 
with that of the SFE; figure 10b shows the results. Here, 
we find a remarkable similarity in the radial profiles of 
i?3_ 2 /i-o and SFEs. Both of them show a strong peak 
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Fig. 11. The plot of the -R3_2/i_o vs. SFE in each annuli 
of M 83. A good correlation can be seen between i?3_ 2/1— 
and SFEs. 

at the center, whereas no significant enhancement can be 
seen at the bar-end, despite the fact that there is a sec- 
ondary peak in the CO( J =1 — 0) intensities, correspond- 
ing to total mass of molecular gas, and the SFRs (figure 
7a and figure 10a). Further, figure 11 shows a correlation 
between R3-2/1-0 an d the SFEs. These suggest that the 
physical state of molecular gas is intimately linked to the 
property of massive star formation. 

Let us compare the properties of star formation at the 
center and those at the bar-end. At the bar-end, the 
amounts of molecular gas and massive stars are enhanced 
compared with other disk region; however, there is no ex- 
cess of SFE. This means that the enhancement in Esfr 
at the bar-end is solely due to an increase in the amount 
of matter (i.e., molecular gas). In contrast, we see a strik- 
ing enhancement in the efficiency of star formation, in 
addition to the elevated molecular gas mass and the rate 
of massive star formation in the nuclear starburst region. 
This means that the presence of nuclear starburst (i.e., 
the striking excess of star formation) is not only due to 
the enhancement in gas mass but also due to the enhance- 
ment in the efficiency of star formation. In other words, 
a simple summation of star forming regions at the bar- 
end (and the disk region) cannot reproduce star forming 
regions in the nuclear region of M 83. 

What causes this difference between the star formation 
properties in the nuclear starburst and disk star forma- 
tion in M 83? It is intriguing to see evidence that both 
i?3_2/i-o and SFE follow the same trend (figure 10b and 
figure 11). This may lead us to speculate that the spatial 
variation of a dense gas fraction governs the spatial varia- 
tion of SFE in M 83; there might be a threshold fraction, 
and once a dense gas fraction exceed the critical value, the 
SFE then increase there. 



Of course, there is an opposite possibility; an excess of 
i?3_2/i-o could be a consequence, not a cause, of the ele- 
vated efficiency of star formation at the center. The strong 
UV radiation and successive supernovae within a starburst 
region could drastically change the physical properties of 
the molecular gas. Further mufti- line/transition observa- 
tions and their modeling using LVG and PDR codes will 
be required to understand the origin and interpretation of 
the observed #3-2/1-0 an d SFE correlation in M 83. 

5. Summary 

We have performed the CO( J = 3 — 2) emission obser- 
vations with the ASTE toward the nearby barred spiral 
galaxy M 83. The size of CO(J = 3 — 2) map is about 5' 
x 5' (6.6 x 6.6 kpc), including the entire inner disk. In 
this paper, we focused on the radial profiles of molecular 
gas properties and their relation to star formation. The 
summary of this work is as follows. 

1. Significant emission was detected at the nucleus, 
bar, and inner spiral arms. We successfully resolved 
these structures at a resolution of 22". No signif- 
icant emission lines were detected in the inter-arm 
region. Our CO (J = 3 — 2) maps spatially correlated 
well with 6 cm radio continuum map, corresponding 
to the distribution of recent star formation. 

2. The luminosity of CO( J = 3-2) I/ CO ( 3 _ 2 ) is 9.3 x 
10 7 K km s -1 pc 2 at the center (r < 0.5 kpc) and 
4.2 x 10 8 K km s _1 pc 2 in the disk region (0.5 < r 
< 3.5 kpc), respectively. The global ^co(3-2) was 
5.1 x 10 8 K km s _1 pc 2 . This means that £co(3-2) 
in the disk region significantly contributes to the 
global I/ co(3 _ 2) . 

3. The average -R3-2/1-0 a t the 22" resolution was ~ 1 
at the center of M 83 (r < 0.5 kpc). The ratio drops 
to a constant value, 0.6-0.7, through the disk region 
(0.5 < r < 3.5 kpc). 

4. The radial profile of the SFE, determined from 6 
cm radio continuum and the CO (J =1—0) emis- 
sion, shows the same trend as R3-2/1-0] i- e -> the 
SFE shows a strong peak at the nucleus (r < 0.25 
kpc), whereas it drops to a constant value in the 
disk region (0.5 < r < 3.5 kpc). 

5. At the bar-end of M 83 (r ~ 2.4 kpc), the amounts of 
molecular gas and the massive stars are enhanced, 
whereas there is no excess of #3-2/1-0 an d SFE in 
that region. This means that the presence of nuclear 
starburst is not only due to the enhancement in the 
gas mass but also due to the enhancement in the ef- 
ficiency of star formation. In other words, a simple 
summation of the star forming regions at the bar- 
end (and the disk region) cannot reproduce the nu- 
clear starburst of M 83. These results could suggest 
that the spatial variation of the dense gas fraction 
traced by i? 3 _ 2 /i_o governs the spatial variation of 
SFE in M 83. 
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